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Abstract 

We study a network model of two conductance-based pacemaker neurons of differing natural 
frequency, coupled with either mutual excitation or inhibition, and receiving shared random in- 
hibitory synaptic input. The networks may phase-lock spike-to-spike for strong mutual coupling. 
But the shared input can desynchronize the locked spike-pairs by selectively eliminating the lagging 
spike or modulating its timing with respect to the leading spike depending on their separation time 
window. Such loss of synchrony is also found in a large network of sparsely coupled heterogeneous 
spiking neurons receiving shared input. 

PACS numbers: 87.19.1m, 05.45.Xt, 87.19.11 
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Phase locking or synchrony of oscillatory systems is a ubiquitous phenomenon in physics 



, l2| , chemistry [Sj , and biology 



m 



The theory of weakly coupled oscillators predicts that 



phase-locking can arise from mutual interactions among oscillators 
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Phase locking is also facilitated by a common external input |14j . 1151 ] either with or without 
weak mutual interactions. But the theory of weakly coupled oscillators is strictly applicable 
only for oscillators with similar natural frequencies 16j. In some brain structures, there are 
mutually-interconnected oscillatory neurons with more frequency heterogeneity than allowed 
by the theory of weakly coupled oscillators [17]. Oscillators with substantially different 



natural frequencies may still become phase-locked if coupled strongly enough 18 



but the 



influence of common inputs to these strongly coupled systems has not been fully explored 
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221 ]. We are specifically motivated by studies of the subthalamic nucleus (STN) 



and globus pallidus (GP) which contain interconnecting oscillatory neurons that receive 
common inhibitory inputs. These neurons are normally uncorrelated but become strongly 
correlated in Parkinson's disease [3]. In this Letter we construct oscillating two-neuron 
network models with different natural frequencies but that are phase-locked due to strong 
mutual interconnections. We show that a common external inhibitory input can reduce or 
abolish phase-locking and decorrelate their firing. The amount of decorrelation depends on 
the phase difference between the spike times in the phase-locked state. 

Our results are demonstrated with coupled two-neuron networks that use three different 
models of varying complexity: a detailed neuron model that reflects realistic properties of 
the STN neurons including rebound currents, a reduced three- variable model (N-l-K) derived 
from the first model, and the generic Hodgkin Huxley (HH) neuron model, which was also 
used as a model of the GP neuron. The results are also extended to a larger network of 
10 sparsely coupled neurons. We used a shared external inhibitory input, because both 
the GP and STN receive strong inhibitory inputs. Interconnections among neurons within 
the STN network model are excitatory, whereas the GP cells are coupled by inhibitory 
connections. Unlike integrate-and-fire type neuron models that may fail to synchronize for 
strong coupling [24,] or for common input [l^, our neuron pairs can display robust 1:1 phase- 
locking of spike times in the absence of external input. This kind of 1:1 phase- locking is 
sometimes also called linear synchrony {25! . [2^ in studies of coupled chaotic oscillators. Our 
study highlights the importance of the timing of the arriving inputs in decorrelating the 
phase-locked state. However, we confine our study to neuron models that are periodic and 
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are not chaotic. 

We consider two coupled pacemaking model neurons at dissimilar frequencies {001^2) and 
affected by identical synaptic noise. The current balance equations (with unit membrane 
capacitance) are 

V, = h-Il- gS2 X {V, -E)- I,^^{Vu t), 

V2 = h-Il- gSl X {V2 -E)- /i,h(V2, t), 

where V1.2 are the membrane voltages. For mutual excitation, E = E^xc = 10 mV and 
g = (yiexc, and for mutual inhibition E = E'inh = —85 mV and g = gi^^. Ij is a steady applied 
current, j = 1, 2 is the index of the neuron. All the three models were used with mutual 
excitation, and the HH model was also used with mutual inhibition. For the realistic STN 
model, = IgiVjjUj, hj,rj, [Ca]^), is a function of Vj, potassium activation (rij), sodium 
inactivation (hj), low-threshold calcium inactivation {rj), and calcium concentration ([Ca]j) 
variables. These two equations along with those for rij, hj, rj, and [Ca]j, together describe 
the coupled system [sG]. For the HH model, = Ij{Vj, rrij, hj, rij), and the gating variables 
are given by standard equations [27] at T = 0°C, Qio = 3('^-6.3)/io_ transformation 
V —V — 60 was used, and E\^ was depolarized from —49.387 mV to —17 mV to trigger 
pacemaking. The N+K model is derived from the STN by setting rebound (Gt) and all other 
calcium dependent maximal conductances (Gca, Gahp) to zero, and thus = /^(V^-, n^, hj). 
The conductance state variable, 5*1,2 = Siyx^A) is noninstantaneous and depends on the 
time history of Vi,2 via a differential equation |28|, a nd can range from 0, when the membrane 
is near rest, to 1 during an action potential 371]. This incorporates a mutual coupling 
time constant, /inh(V^',^) is an inhibitory input current generated from alpha function 
conductances timed at random but identical Poisson process arrival times (tj, i = 1, . . .) with 
rate Iinh{Vj, t) = dnh Eu 7^^^ '-'^ H{t - t^) x {Vj - E-,^^). Gi„h = 0-1 nS/zim^ for 

STN and N+K, and — 1 mS/cm^ for the HH network, and Tjnh =1 — 10 ms. A fourth order 
Runge-Kutta algorithm with 1 /xs time step was used to integrate up to 1000 s. The spike 
time sequences (si and S2) of our networks, when they were not driven, are 1:1 phase-locked 
with spike separation and thus fall in the (si, S2) phase space on the diagonal line si = 82- 
The deviations from this linear synchrony (also see [251 and references there in) due to shared 
input are large and random. An asynchrony parameter that measures the fraction of such 
deviating spike times is defined: A = 1 — A^sync/^max, where Agync is the number of spike time 
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FIG. 1: Emergence of asynchrony for shared inhibition, (a) Coupled STN network, toi = 2.7 Hz 
(/i = 0) and UJ2 = 9.3 Hz {h = 8 pA/fim^). dnh = 1 nS/^m^. ^r*^^ = 0.4 nS/^m^. The top inset 
shows brief time courses at Qexc = 0.5 nS//im^. For the same conductance parameters, the effect 
of frequency heterogeneity is shown in the bottom inset, (b) Frequency dependence of A for the 
three networks: STN (Gjnh = 1 nS//um^, and three different values of ^exc), N+K (/i = 0, /2 = 8 
pA/^m^, 5rexc = 0.62 nS/fim^, Ginh = 1 nS/^m^), and HH {h = 0, h = 8 /iA/cm^, ^exc = 0.2 
mS/cm^, Ginh = 1 mS/cm'^). pi, . . . , mark input frequencies at which A changed its character 
in the STN network due to T-current activation, and are described in the text, (rinh = 1 ms). 

pairs (one spike from each neuron) in a time window t^^ and A^max = maxjA^^i, A^2}, where 
are spike counts of the two neurons. A = for 1:1 phase-locked state. For phase-locked 
states that have large [in Fig. EJ^b)], shorter separations also constitute deviations from 
;he coherent state. Thus an asynchrony measure called inter spike interval (ISI) distance 
29] defined as j; Jq J'^^^sl^l^y dt is used. The ISI-distance quantifies the jitter around the 
diagonal line, si = S2. However, both these measures give qualitatively similar results. 

Figure [It^a) shows the results of two mutually excitatory STN neurons that have very 
different intrinsic spike frequencies. In the absence of external input, mutual coupling that 
exceeds a critical level ((^cxc) phase-locks the spike times of both the neurons in 1:1 ratio. This 
causes A to become as the frequencies of both the cells become identical. An example 
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of the voltage time courses in this state is shown in Fig. [TJ^a), top inset, left. When the 
external inhibitory input was turned on [Fig. [T](a), top inset, right], the spike times of both 
the cells became arrhythmic and the firing rates became different. The common input 
countered the effect of mutual coupling by decreasing the number of the 1:1 phase-locked 
spike pairs. The enhanced spike count resulted due to activation of low-threshold rebound 
calcium current (i.e. the T-current). But the loss of phase- locking is not dependent on the 
T-current activation (see results using the N-l-K network below). The result of the loss of the 
1:1 phase- locking between the spike times is also seen as a finite non-zero value of A in the 
earher phase-locked regime {g > (^exc)- this network, frequency heterogeneity is essential 
to elicit finite A. This is demonstrated in Fig. [I](a), bottom inset. In the absence of input, 
phase-locking was achieved for a wide range of frequency disparity. With inhibitory input, 
A became non-zero in this region except when the frequencies were identical (Ji = 12). The 
dependence of A on the input rate is shown in Fig. W(Jo) for all the two-neuron coupled model 
networks with mutual excitation. In the absence of input (i.e. at Amh = 0), the networks 
were 1:1 phase- locked (i.e. A = 0). At an input rate of 1 kHz, about 40% of the locked 
spikes in the N-l-K network, and more than 50% of the locked spikes in the HH network 
became unpaired. But the STN network model exhibited frequency-dependent activation of 
its rebound current that led to a modulation of A. The T-current is most actively recruited 
in the slower neuron between pi and p^, and in the faster neuron between p2 and ^4. Thus, 
the firing rate as a function of Amh increased in both the neurons between p2 and ps leading 
to an increasing A. Active recruitment in only one of the cells led to enhancement of phase- 
locking or reduction in A. 

The loss of 1:1 phase- locking observed in the N-l-K, the HH, and the STN network (for 
Ainh < Pi, hence no significant rebound spikes) emerges by selective elimination of one of the 
spikes in otherwise phase-locked spike pairs [see ISI histogram in Fig. EJ^b) bottom inset]. 
Shared input cannot desynchronize the locked state if the spike pair separation = 0. For 
small tw, a brief synaptic input is needed for disruption of the locked state, tw is determined 
by the dynamics of the synaptic coupling conductance (S'(V^-,t)) and the coupling strength 
(7exc- It increases with frequency heterogeneity [Fig.[2](a)]. For the STN network, empirically, 
tw ~ CL{gcxc — dcxJ^^j where a = 0.78, b = 0.44 away from the criticality {g*^^) ? ^ — 1-3) 
b = 0.2 near criticality. An inhibitory stimulus event preceding a locked pair [Fig. [2](a) inset] 
may reduce the spike separation, but one or more events that are positioned appropriately 
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FIG. 2: (a) Time window of spike separation in 1:1 phase-locked state in the absence of extrinsic 
input for the coupled STN, N+K and HH excitatory neuronal pairs as a function of frequency 
heterogeneity (Ii is fixed at 4 pA/^um^ for STN and N+K, and at 4 /^A/cm^ for the HH.) The 
inset shows scenarios of input events with respect to spike pair separation discussed in the text, 
(b) Excitatory STN network. Predicted and computed time windows (see text). The top inset 
shows the mechanism by which oncoming inhibition suppressed a lagging spike. The bottom inset 
shows the resultant interspike interval histogram densities, (rjnh = 1 ms, ^exc = 0.5 nS/fira^.) 

in a time window tp (smaller or bigger than t^) with respect to, but preceding the onset of 
the excitatory postsynaptic conductance, can eliminate the stimulated spike in the lagging 
neuron. A is a measure of the probability of such a favorable window receiving one or 
more inhibitory stimuli (1 — Po('^inh, ^p), where P„(A,t) is the Poisson probability density) 
that successfully eliminate the following spike. Assuming a successful elimination of the 
secondary spike by one or more arriving inputs in a window, and assuming tp ^ t^, we write 
A = Q{1 — e"'^'"'^*"'), where Q is a constant that depends on the phase- locked frequency. 
(iy4/(iAinh = QtuiC"^'"''*'". Q was estimated numerically from dA/dXi^h- At Gjnh = 0.1 
nS//im^ and = 1.8 ms, Q = 0.034. By using Q, and the numerically evaluated dA/dXi^i,, 
we solved the above transcendental equation for at several Ginh values, and plotted it as 
tp in Fig. [2](b), in dots . A similar window for suppression of the lagging spike is computed 
deterministically by using an inhibitory conductance input preceding the phase-locked spike 
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FIG. 3: (a) Dependence of A on Tinh- For STN network Ainh = 100 Hz, g^xc = 0.5 nS/fiia^. For 
HH network Gjnh = 0.5 mS/cm^ and gexc = 0.1 mS/cm^. Ii and /2 are as in Fig. [TJb). 1//? = 3 
ms for both the networks, (b) Mutually inhibitory network of two coupled HH neurons (/i = 0, 
/2 = 8 ^A/cm^) displaying spike time asynchrony for inhibitory input. At Ainh = 0, t^, = 7.8 for 
#1 and 9.6 ms for ^2. Voltage traces transitioning from phase-locked (Ainh = 0) to asynchronous 
state for finite input rate are illustrated in the top inset. The phase response curves corresponding 
to the two phase-locked neurons are illustrated in the left inset, bottom. The right inset shows the 
densities of spike pair separation in the asynchronous state. 

pair, and is shown in Fig. Mjo) (thick solid line) for comparison. Alternatively, A may be 
approximated. Such an approximation even at large Ginh values well predicts the growth 
of A [Fig. [U^b), STN, analytic] at small Ainh: A = 1.02 (1 - e-ooo^Vh), Ugjng a similar 
approximation for the HH network [Fig. [T](b), HH, analytic] yields A = 1.8 (1 — e^° °°^^'^'"h^_ 
For basal ganglia networks, values of Tjnh are near 10 ms. Increasing Tinh generally decreased 
the firing rates in the models because the input turned increasingly deterministic, and the 
asynchrony parameter showed a moderate increment [Fig. [3]^a)]. But the STN network 
showed up to 80% loss of phase-locked spike pairs when Ginh is reduced appropriately at 
long Tinh such that the cells could maintain high firing rates due to activated T-current. The 
modulation in A is similar in nature to that seen in Fig. [T](b). 

We now briefly discuss the behavior of mutually inhibitory HH neurons [Fig. [3](b)] in 
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response to common external Poisson inhibition. Phase-locking can now occur in such 
networks with near the anti-phase state, is modified by the frequency heterogeneity. 
This phase difference between the spikes is either decreased or increased by the occurrence of 
oncoming inhibitory synaptic conductance events before or in between the phase-locked spike 
pair. This increased the asynchrony as measured by the ISI-distance (defined earlier) with 
increasing Amh at fast as well as slow synaptic time constants, Finally, we illustrate the 
manifestation of these two forms of asynchrony in a larger network (A^ = 10) of excitatory 
[STN, all-to-all coupled; n = 9; Fig. |4](a)] and inhibitory [HH, sparse coupled; n = 4; 
Fig. 11(b)] neurons. The current balance equations are Vj = Ij — /J — Y^=\ MjkSk x (Vj — 
E) — Iinh(yj,t), j = 1, . . . , 10. Mjj = 0. Mjk = 1, if the presynaptic neuron k is connected 
to the postsynaptic neuron j, and otherwise. In both the networks, in the absence of the 
input, the spike times were asynchronous for weak coupling (see ^o| for T-current induced 
clustering), but strong mutual coupling above a critical level synchronized the average firing 
rates of all the neurons to a common frequency. When the shared inhibition is turned on, 
the STN network neurons fired faster due to T-current activation, and the HH network 
neurons fired slower. But the average firing rates exhibited asynchrony by diverging from 
the common locked frequency. 

In contrast to the earlier studies on weakly coupled oscillator theory, we showed that in 
strongly coupled heterogeneous networks, shared inhibition can lead to asynchrony of spike 
times. Our simulations also revealed similar role for shared excitation. Our results lead 
to the possibility that modulating the synaptic input strength or input frequency might 
lead to switches in the synchrony-asynchrony (or coherence-incoherence) transitions in cou- 
pled oscillatory networks that receive coherent input. It would be of interest to study the 
emergence of asynchrony due to shared inputs when the phase-locked state is affected by 
background noise 31| or if the oscillators themselves behave chaotically 32|- In Parkinso- 
nian network models of deep brain stimulation, earlier studies 33|, l34] relied on the network 
connectivity between different brain nuclei to produce complex asynchronous patterns. Our 
results suggest of an alternative mechanism for generating such asynchrony within each brain 
nucleus. 

We thank the Computational Biology Initiative (UTHSCSA/UTSA) and Texas Advanced 
Computing Center (TACC), The University of Texas at Austin for providing high perfor- 
mance computing resources. The work was supported by NIH/NINDS NS47085. 



8 




FIG. 4: Emergence of asynchrony from coherent states in networks of = 10 excitatory all-to-all 
coupled STN (a) and sparsely coupled inhibitory HH (b) neurons measured by average spike rate 
as a function of mutual coupling without and with shared inhibition, (a) Ij = {j — 1)10/9 pA//Ltm^, 
j = 1, . . . ,N, and Ginh = 0.5 nS//xm^. (b) Inset shows My (with n = 4); a filled square represents 
a synaptic connection. Ij = 5 + {j — l)8/9/xA/cm^, j = 1, . . . ,N, Ginh = 0.1 mS/cm^. 
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